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Abstract. Within the framework of the LSW theory of Ostwald ripen- 
ing the amount of volume of the second (solid) phase in a liquid solution 
that is newly formed by recrystallization is investigated. It is shown, 
that in the late stage, the portion of the newly generated volume formed 
within an interval from time to to t is a certain function of t/to and an 
explicit expression of this volume is given. To achieve this, we introduce 
the notion of the return radius r(t,to), which is the unique radius of a 
particle at time to such that this particle has - after growing and shrink- 
ing - the same radius at time t. We derive a formula for the return radius 
which later on is used to obtain the newly formed volume. Moreover, 
formulas for the growth rate of the return radius and the recrystallized 
material at time to are derived. 



1. Introduction 

Recrystallization of minerals is a combination of simultaneous processes 
of dissolution and precipitation that leads to formation of larger mineral 
crystals. A driving force of recrystallization in geological environments is 
usually either the difference between lithostatic and hydrostatic pressures or 
the dependence of the chemical potential of the interface on the grain size. 
This type of coarsening process, where larger particles are growing at the 
expense of smaller ones is called Ostwald ripening. A prominent example is 
Ostwald ripening of calcite (CaCOs) in an aqueous solution, which is subject 
of a number of experimental investigations, see e.g. [5], [7], [3],[10],[2],[1]. 

In these experiments one is interested in the amount of newly formed 
crystalline material, i.e. the volume of the solid phase present at time t, 
that has been precipitated from the solution after some time instant to < t- 
Experimental methods tending to determine the newly formed material are 
usually proceeding as follows: An isotope of a crystal lattice constituent 
of the mineral - in case of calcite Ca or 14 C - is added as a radioactive 
tracer to the solution at time to- The isotope concentration in the solution is 
monitored during the experiment, and changes of the isotope concentration 
happen due to two processes: either it is a diffusion into the crystal, or 
burial in the calcite layers newly formed during recrystallization. 

In order to derive the amount of newly formed material from the Iso- 
tope uptake, in [1],[3],[2] it was argued that the amount of isotope uptake 
is proportional to the newly formed material. However, this approach does 
not account for the surface changes during ripening. Whenever the surface 
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area decreases, adsorbed isotope will be returned to the solution which does 
influence the overall isotope uptake. The question we were asked by exper- 
imentalists was if, and under which conditions, the surface change may be 
neglected or, more precisely, how the influence of surface change compares 
with the amount of newly formed material concerning the isotope uptake. 
To answer these questions, a first step is to derive analytical expressions 
both, for the newly formed volume and for the change of the surface area. 
The former will be the subject of this paper, while the latter usually may 
be derived very easily in some mean field theory. We will not attempt to 
draw conclusions concerning experiments, since this would require too many 
details of geochemistry in order to construct a sound mathematical model. 
This is planned to be done in a forthcoming article. Besides these moti- 
vations from experiments, an explicit formula for the newly formed volume 
during Ostwald Ripening is also interesting from a theoretical point of view 
in order to better understand the long time behavior of the coarsening sys- 
tem. 

In this paper we are going to calculate the newly formed volume theoreti- 
cally in the late stage of the coarsening process, assuming that the so called 
LSW theory is valid, see [6], [9]. In this mean field theory, it is assumed, 
that the grains are spherical particles and the growth kinetics of a grain only 
depends on its size compared to the size distribution of all particles and not 
on the local environment. Moreover it is assumed, that the total volume of 
the crystallized material is conserved, i.e. large particles grow exclusively at 
the expense of smaller ones.. A short review of the LSW theory is given in 
the next section. In Figure 1 we have sketched how the size of particles may 
have changed from time to to t. Note that a particle may grow at time to 
but start to shrink at some t\ with to < t\ < t leading to R(t) < R(to), see 
Figure 2, where some typical examples of the time evolution R(t) of particle 
radii R, are shown. 



Figure 1. Newly formed solid phase. During Ostwald ripening, 
bigger particles are growing at the expense of smaller ones. Com- 
paring the sizes at time to and at time t > to, the three pos- 
sible cases are sketched: Ri(t) > Ri(t ), i? 2 (*) = ^2(^0) an d 
Rsit) < Ra(to)' Thus at time t > to, a certain amount of the 
solid phase, marked in gray, has been newly formed between time 
to and time t. 
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t(r)/t = 5 




t/t ,t Q = 1000 



Figure 2. Time evolution of the radius R(t) of a spherical particle 
with different initial conditions Rq — i?(io) in the late stage of 
Ostwald ripening in the diffusion limited LSW regime as described 
in section 2, with non dimensional critical radius i? c (0) = 1 . Since 
in all three cases the initial radius R(to) is larger than the critical 
radius R c (to), the particles start to grow at the expense of smaller 
particles. However, after a certain time they start shrinking at the 
expense of even bigger particles, since the critical radius is also 
growing in time. For a particle radius i?(io) — being larger 
than the critical radius at time to, there is a return time t(Ro) 
such that R(t) = Rq. The other way round, given t > to, there is 
a unique radius r such that a particle with initial radius R(to) = r 
will have the same radius at time t, i.e. R{t) = r. This radius 
will be called the return radius. E.g. for the blue line we observe 
that i?(2<o) = R(to), i.e., the return time is t = 2to and the return 
radius for t = 2to is given by R(to) rj 8.42. In fact, we have 
calculated suitable initial values R(to) for the above cases with our 
method described in section 3. 



For an ensemble of particles with radii Ri(to), the newly formed volume 
V new (Mo) between t and t is 

(1) y ncw (Mo) = h E (^W 3 " ^(*o) 3 )- 

i:R % {t)>Ri(t ) 

Note that due to the assumption of mass conservation of the second phase, 
the newly formed crystalline material all comes from dissolution of smaller 
particles. In a mean field picture Ri(t) > Rj(t) if and only if i?j(to) > 
Rj(to). As will be discussed in more details in the next sections, there is a 
critical radius R c (t) such that at time t precisely the particles with radius 
R(t) > R c (t) are growing. However, R c (t) is growing faster then R(t) and 
particles that initially grow, start to shrink at later times, see Figure 2. 
Thus, there is a unique radius r = r(t,to), such that i?j(to) = Ri(t) if and 
only if Ri(to) = r. We will call this radius r the return radius and V new (t, to) 
may be expressed as 

(2) v ncw (Mo) = h E (^(*) 3 - R ^f)- 

i:Ri(t)>r(t,t ) 
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We will show in section 3 that in the late stage of Ostwald ripening in the 
LSW regime the return radius r(t, to) is a function of t/to and can be easily 
computed by inverting an explicitly given function. Moreover we obtain an 
analytic expression for the growth rate of r(t, to) at t = to. This will lead 
in section 4 to our main result: Also the volume V new (t,to) or equivalently 
the volume fraction <3? new (t, to) of the newly formed solid phase within the 
interval from to to t depends on t/to and may be calculated by function 
inversion of explicitly given expressions. We also give an explicit expression 
for the initial rate of the formation of new solid material. 



We will shortly review the main results of the LSW-analysis given in 
[6], [9], for a more detailed description see e.g. [8]. In the following the three 
dimensional model is presented, using the notations of [4], where the two 
dimensional case is discussed. The kinetics of Ostwald Ripening is governed 
by two different processes: the mass transport between the mineral grains 
via diffusion in the solution and the attachment/detachment process at the 
grain boundaries. The two limiting kinetic regimes are termed attachment 
limited (AL) growth - here diffusion is assumed to be infinitely fast compared 
to the attachment process - and diffusion limited (DL) growth - assuming 
the attachment process to be instantaneous. 

The first assumption of LSW-theory is, that the volume fraction of the 
dispersed solid phase is very small. Thus it is assumed that the solid phase 
consists of many disconnected particles far away from each other which 
moreover are assumed to be spherical and have immobile center. For a 
large system, this ensemble of particles may be characterized in terms of a 
particle radius distribution function F(R,t). The number of particles per 
unit volume is then given by n(t) = J Q F(R, t) dR. Assuming that no 
nucleation and coalescence of particles takes place, F obeys the continuity 
equation 



Note that the number of particles n(t) may only change due to radii de- 
creasing to zero in finite time. 

The second assumption is, that far away from the particles, the chemical 
potential u(x,t), describing the change of energy per change of mass of the 
dissolved material may be approximated by a spatially constant mean field 
u(t). Finally, one assumes that the volume fraction $ is constant in time. 
The driving force of Ostwald ripening is the difference between the chem- 
ical potential u(t) in the solution and the chemical u e q(R) at the particle 
surface, depending on the curvature 1/R. From the above assumptions it 
is derived, that the growth rate R(t) of any particle may be described in 
non-dimensional form as 



2. LSW Theory for Ostwald Ripening 



(3) 



d t F + d R (RF) =0. 



R(t) 



( 



u(t) 



U eq (R)) 



1 



u(t) 




(DL) 



R 



R(t) 



(u(t) 



u eq (R)) 





(AL). 



E^ 2 
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Moreover, assuming the chemical potential at the surface to be given in non 
dimensional form by the Gibbs Thomson law u cq (R) = 1/R, the critical 
radius R c such that particles with radius R < R c are shrinking while the 
ones with R > R c are growing is given by 

R c (t) = l/u{t) 

and the growth law may be written as 

1/R ,\ 1 



(4) m = W { K ~l), fl(*)=j (DL) 

(5) R(t) = U^-l), fl(t) = 4 (AL). 



r\r c r w R2 

Introducing the new variables 

(6) Z =^ T = ln W' 

eqs. (4), (5) become 

dz z — 1 

(7) 



dr z x 



where A = 2 (DL) or A = 1 (AL) and u is a function of the critical radius 
R c : 

(8) v= ! , (DL) v= 1 , (AL). 

R 2 c R c (t) R c Rc{t) 

Note that in the new coordinates, z = 1 corresponds to the critical radius 
R c (t) of a particle, that is neither growing nor shrinking. For later use we 
also point out, that according to (7) the rescaled radius z may be shrinking 
also for values z > 1, i.e. for growing particles. Actually, as will be seen 
later, the growth rate dz/dr is negative for all particles in the late stage. The 
essential point of the LSW-analysis is, to argue, that v becomes constant at 
late times approaching the unique values ([9], [6]) 

(9) ^ = 27/4 (DL), v = \ (AL). 

This in turn implies a scaling law for the critical radius R c by solving the 
differential equations for R c (t) given in Eq. (8). Assuming that t = is 
already in the late time regime, one does obtain 



R c {t) 



{Rc(0) 3 + lt) 1/3 (DL) 
(i? c (0) 2 + it) 1/2 (AL). 



For later purpose we note that this scaling law implies for i, to S> i? c (0) 3 , 
that 

(10) ^M = ( t ^Y / \ 7 = 3(DL), or 7 = 2(AL). 



R c (t) V t 
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Moreover, z(t) may be obtained from (7) by inverting the explicit solution 
of t(z) given by 

- 1 \-l 



r(z) 



(11) 




dz 



- § ln(z + 3) - § ln(3 - z) 
ln(2 - z) 



(DL) 
(AL)' 



where we have omitted a constant of integration, since later on only differ- 
ences t(z) — t(zq) will be used. Here the constant values of v given in (9) for 
late times have been used. In Figure 3 the solution z(t) with initial value 
z(0) = 5/4 is depicted. Note that the reduced radius z(t) is decreasing in 
time also for particles with z > 1, i.e. R(t) > R c (t). In fact, one easily veri- 
fies, that the right hand side of (7) is negative for all particles with z < 3/2 
(DL) or z < 2(AL), respectively. 




Figure 3. The rescaled radius z(t) for a particle with initial value 
z(0) = 5/4 obtained by inverting t(z) given in (11). Note that de- 
spite the fact, that the initial particle radius is larger than the crit- 
ical rescaled radius z = 1, the rescaled radius z(r(t)) = R(t)/R c {t) 
is decreasing in time. 



(12) 



The continuity equation (3) for f(z,r) := F\R c z,t(T)\ R c , 

dz 



0, 



can now be solved by a separation ansatz f(z,r) = g(r)h(z), yielding the 
scaled normalized island size distribution functions, see [9], [6] 



(13) h(z) 



'81e2-§z 2 (z + 3)-5(§ - z)"ir exp (gf^) : z< 
: z > 



(14) h(z) -- 
and the scaling 
(15) 



24z(2 - z)- 5 exp (f^f) : z<2 
: z > 2 



9{t) = 9oexp(-3r) = g 



(DL) 
(AL) 
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Note that one outcome of the LSW analysis is a cutoff of the particle size 
distribution function with the largest rescaled radius given by 

(16) z max = 3/2(DL), z max = 2(AL). 

Finally we remark that the critical radius R c may be expressed in terms 
of the mean radius R as R c = R and R c = |fi for the (DL) and the (AL) 
ripening, respectively [9]. 



3. Return radius 

As has been sketched in the introduction, to calculate the newly formed 
volume between time to and time t we need to calculate what we have called 
the return radius r = r(t,to), i.e. the unique radius r such that a particle 
with radius R(to) = r will have the same radius at later time t, R(t) = r. 
In fact, this amounts to solve a boundary value problem for the differential 
equations (4), (5). We will see, however, that the return radius may be 
calculated quite easily in the asymptotic regime of LSW, i.e. for to large 
enough, without explicitly solving the boundary value problem. Since in the 
following to is fixed, we will write r(t) instead of r(t, to). 

We will use the rescaled coordinates z,t as introduced in (6). First note 
that in the rescaled coordinate z the return radius r does take two different 
values at time to an d at time t, namely z(to,r) and z(t, r). We will use the 
notation To = r(to), r = r(t), zq = z[tq) = z(to,r), z = z{t) = z(t,r). From 
(6) follows, that r being the return radius is equivalent to 



(17) 



In z — In zq = tq — T. 



Moreover, the scaling law (10) implies for large to 

z(t,r) R c (t ) 



(18) 



z 

ZQ 



z(t Q ,r) R c (t) 



The two identities (17) and (18) uniquely fix the return radius r = r(t,to), 
which may be calculated as follows. Since the time dependence of the 
rescaled radius z = z(t) is given by equation (7), we use the explicit so- 
lution of t(z) given in (11) to express (17) as 



(19) 



a(z) = ot(zo), with a{x) = hxx + t(x). 



Since the return radius has to be larger than the rescaled radius, r > R c (to), 
we may assume zq = z (to,r) > 1. One easily checks that 



a'(z) < 



> 
= 
< 



and 



ze (0,1) 
z = l 

Z £ ( 1 j ^max 



lim a(z) 

z—>0 



lim a(z) = —oo. 



Thus, equation (19) allows for a unique solution for zq £ [1, z mgx ). Denoting 
this inverse of a on (— oo, a(l)] by tp, we may express z = z(t, r) as a function 
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(20) 



P ■ [Mmax] -> [0,1] 

p(z ) := ip(a(z )), z £ (l,z max ) 
p(l) = 1, p(z max ) = 0. 



This function is easily evaluated numerically, e.g. using bisection, see Fig- 
ure 4. 





1.2 



1.4 



1.6 



1.8 



-DL- 
-AL 



Figure 4. (Left) The function a(z) = lnz + t(z) in (19) has a 
unique maximum at z = 1. (Right) Choosing z > 1, there is a 
unique z := p(zq) < 1 with a(z) = a(zo). The function p has been 
evaluated numerically using bisection. 



Now we may calculate the return time t(r) for any initial radius r, i.e. 
the time t, such that R(t) = R(to) = r. By equation (18) we have 

(21) 4 ( r ) = t (_fo_) 7 z = r/R c (t ). 

Inverting t = t(r) in (21) numerically yields the return radius r(t), see 
Figure 5. We point out, that according to (21), the return radius r(i,io) is 
a function of t/to. As it will turn out in the next section, we only need the 
pair zq,p(zo) in order to calculate the amount of material which has been 
produced by recrystallization, see Figure 6. Let us finally calculate the 
growth rate r(t) of the return radius at time t = to- Using (21), we get 

'M^M)" 1 = (;£)" Vo). 

and equation (21) gives 

( 22 ) = £ / z Q \7-i p (zq) - z p'(z ) 



dz Q \p{zo)J p(z ) 2 

To proceed further, the derivative of p(zo) defined in (20) is needed at zq = 1. 
Since a'(l) = 0, a"(l) / 0, we have 

(23) p(l + h) = 1 -h + C(h 2 ), i.e. p'(l) = -1. 

Evaluating (22) at zq = 1, and using p(l) = 1 yields 

dt 

24 — = 2 7 t . 

dz o , =i 
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Figure 5. (Left) Return time over rescaled initial radius. (Right) 
Rescaled return radius over time. Time is given in units of the 
starting time to. 



1.5 







\ X 
\ *""*■>. 
\ - .» 

s ""■»«. 
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X 

X 

X 

X 

~ ~ ~ 


— z Q =z(t ,r)(DL) 
---p(z ) = z(t,r)(DL) 

— z Q =z(t ,r)(AL) 
---p(z )=z(t,r)(AL) 





l/t o 

Figure 6. The return radius r — r(t,t ) may be most easily cal- 
culated in rescaled coordinates using (21). 



So we end up with the growth rate 
(25) ffo) = 



4. Volume fraction of recrystallized material 
Now let <3? denote the volume fraction of the solid spherical phase, i.e. 

roc 

(26) $(t) = §tt / F(R,t)R 3 dR 

Jo 

and denote by <3? new (t, to) the volume fraction of that part of the solid phase 
at time t that was produced between time to and time t through recrys- 
tallization. Only those particles, which have grown between to an d t, i.e. 
R(t) > R(to), contribute. These particles are precisely the ones for which 
R(t) is larger than the return radius r = r(t,to), or, equivalently, R(to) > r. 
Thus <l? new (t,to) may be calculated as the difference of the volume fraction 
of all particles being larger than r at time t and the volume fraction of these 
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same particles at time to'- 

(27) $ ncw (Mo) = \A I F(R,t)R i dR- I F(R,t )R 3 dR 



Assuming the LSW-theory to be valid, we may express F(R, t) and F(R, to) 
in terms of f(z, r) = g{r)h(z) and f(z, to) = g(ro)h(z), respectively, as given 
in (13) - (15). Since F(R,t) = R c (t)f(z,r) and F(R,t ) = R c (t ) f (z , r ) , a 
change of variables in (27) leads to 
(28) 

rz{t ,r) pz 
<£> llcw (t,t ) = §7r 50 i? c (0) 3 / h(x)x 3 dx = §7r# i? c (0) 3 / h(x)x 3 dx, 

J z(t,r) J p(zn) 



'z(t,r) Jp(zo) 

where the return radius r = r(t,to) may be calculated from equation (21) as 
described in the last section. In fact it is sufficient to solve for zq = z(to,r) 
in order to evaluate (28), i.e. to invert the function t = t{zo) in (21) and 
to calculate z(t, r) = p(zo), see eq. (20). Similarly we may express the total 
volume fraction in the late stage as 

(29) $(to) = $(t) = |vrc/o-R c (0) 3 / h{x)x 3 dx, 

Jo 

Thus the percentage of the volume of the second phase at time t, that has 
been produced by recrystallization between time to and time t, is given by 

(30) 

<p(t,to) := $ new (Mo)/$ = = / h(x)x 3 dx, z 3 := h(x)x 3 dx, 

Z 3 Jp(zo) Jo 

and may be easily computed by determining zo = z(to, r) and p(zo) = z(t, r) 
as described in section 3 and using numerical quadrature. Again we point 
out, that tp(t, to) may be expressed as a function of t/to, since r(t, to) depends 
on t/to only. 

In Figure 7, the specific produced volume ip(t, to) is depicted over the 
normalized time s := t/to- As expected from (30), we observe that ip(s) 
approaches the value (p = 1 for s — > oo. Moreover, the growth of the 
recrystallized volume is nearly linear up to the value (p(s) ~ 1 /4. The growth 
rate of ip(t,to) at time t = to may be calculated using equations (23), (25) 
as follows 



d_ 
It 



<p(t,t ) 



t=t 



d 

dz 



20 



p(.*o) 



h(x)x 3 dx 



d 



zo(t) 



d 



(31) 



= \ h(zo) - h(p(zo))—p(z ) 
h(l) 



ZiY- 



t=t 
1 

-iWo 



ito 



Here in the last equality we have again used that p{\) = 1 and p'(l) = — 1, 
see section 3. For the growth rate, we obtain from (31) and (13), (14) the 
numerical values 



(32) ^(Mo) 



h(l) _ (27e~ 2 2- 8 / 3 /(zHo) « 0.51/t 
t=t ~ Ito^ ~ \l2e~ 3 /(z^to) « 0.62/to 



DL 
AL 



RETURN RADIUS AND VOLUME OF RECRYSTALLIZED MATERIAL 



11 



5. Conclusions and outlook 

Let us summarize our findings: We fix some time instant to 3> 1 i n the 
late stage of Ostwald ripening, where we assume that the LSW theory is 
valid. All particles with radius R(to) being larger than the critical radius 
R c (to), i.e. with renormalized radius 

z = r/R c (t ) G (Mmax] 

start to grow but - except for the particle with maximum radius z max , see 
(16) - will start to shrink at some later time instant. For every such zq we 
can find the return time t(zo), i.e. the time instant at which the particle 
will return to its initial size, from the monotonic formula (see (21)) 

t(zo)=to GS)) 7 - 

The function p(z) is defined in eq. (20) and its graph is given in Figure 4. 
The inverse function Zo(t) is shown in Figure 5. 

A formula for the specific volume ip(t, to) of the newly formed material 
within the time interval [to,t] has been derived for both, AL and DL ripen- 
ing. It is given in terms of the corresponding scaled normalized island size 
distribution function h(z) as 

I fz (t) 

(33) v?(Mo) = / h(x)x 3 dx. 

const Jpfa®) 

The graph of ip is given in Figure 7. Moreover, an explicit formula for the 



1 




Figure 7. Specific volume ip(t, t ) of recrystallized material from 
time to to time t over normalized time t/to in the late stage, where 
LSW theory is valid. <p(t,to) is computed numerically by first 
solving (t/to) 1 / 7 = z o/p( z o) f° r z o = z ( r ,to)i with p given in (20), 
and then calculating the integral S p ? Zo \ h(x)x 3 dx using numerical 
quadrature. The slope at t/t = 1 is given in (32). 

growth rate 4^(t,t )k=t is S iven in ( 32 )- 
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In a forthcoming paper we plan to apply our findings to the isotope uptake 
during Ostwald ripening. Here one also has to account for the adsorption of 
the tracer material at the surface of the solid phase and therefore the change 
of surface area during ripening may play a role. Moreover, a generalization 
to the case of large volume fractions, where the LSW theory is not valid, 
needs further investigations. 
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